scalar DiNum =  tpyrolysis->solidRegionDiffNo();

